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It is shown that in perfectly quasi-isodynamic stellarators, trapped particles with a bounce frequency much 
higher than the frequency of the instability are stabilizing in the electrostatic and collisionless limit. The col¬ 
lisionless trapped-particle instability is therefore stable as well as the ordinary electron-density-gradient-driven 
trapped-electron mode. This result follows from the energy balance of electrostatic instabilities and is thus 
independent of all other details of the magnetic geometry. 


Stellarators seek to confine fusion plasmas by means of a three-dimensionally shaped magnetic field. In recent years, the 
art of optimizing this field to improve plasma performance has taken great strides. In particular, it has proven possible to 
shape the magnetic field in such a way that the collisional (so-called neoclassical) transport is reduced almost to the level of 
axisymmetric devices. An important question that then arises is how this optimization affects the properties of microinstabilities 
and the turbulence they tend to cause. In tokamaks, most of this turbulence is driven by ion and electron-temperature gradient 
(ITG and ETG) modes and by the trapped-electron mode (TEM). In this Letter, we demonstrate that one of the most important 
classes of orbit-optimized stellarators, so-called quasi-isodynamic ones, is automatically immune to the ordinary TEM and to all 
lower-frequency electrostatic instabilities, if the temperature gradients are small enough compared with the density gradient and 
collisions can be ignored. Quasi-isodynamic stellarators could therefore benefit from reduced transport both in the neoclassical 
and turbulent channels. 

A toroidal magnetic field B is quasi-isodynamic when the contours of constant B = |B| are poloidally (but not toroidally) 
closed and all collisionless orbits are perfectly confined lllli . Thus, if B = X Va, where tp denotes the toroidal flux, then 
the radial drift should vanish when averaged over the bounce time Tt, 
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where the integral is taken along the field between two successive bounce points, is then constant on flux surfaces, i.e., J depends 
on Ip, the energy and magnetic moment of the particle but is independent of a. Wendelstein 7-X is the first stellarator to approach 
quasi-iso^namicity, and substantially more quasi-isodynamic configurations have been found computationally in the last few 
years 10-0. 

These devices are so-called maximum- J configurations, where J has a maximum on the magnetic axis and dJ/dip < 0. 
In 1968, Rosenbluth 10 had already noticed that the maximum-J property is beneficial for the stability of interchange modes 
with frequencies above the drift frequency but below the bounce frequency of all plasma constituents (see also Refs. AH). 
Eor typical microinstabilities (except the ETG mode), this condition holds for the electrons but not for the ions. Nevertheless, 
considering the full gyrokinetic system of equations, we show in this Letter that stability prevails far beyond the limit considered 
by Rosenbluth. 

The physical reason for this remarkable stability has to do with the direction of the precessional drift of the trapped particles. 
If the wave vector perpendicular to the magnetic field is = kaVa k^ipVip, we define the magnetic drift frequency ojda = 
kj_ ■ Vda and the drift wave frequency = {Taka/ea)dlnna/dtp for each particle species a in the usual way, where the density 
ria and temperature Ta are constant on flux-surfaces, and \da = b x {{v]_/2)V InS + uj|K:)/f2a denotes the drift velocity, with 

b = B/|B| the unit tangent vector and k = b • Vb the curvature vector of the magnetic field. The precession frequency then 
becomes 


— ka^CX * Vda V tp * \da i 

where an overbar denotes the bounce average. By design the last term vanishes for quasi-isodynamic configurations. The 
remaining term can be expressed as a derivative of the parallel adiabatic invariant, taken at fixed energy and magnetic moment, 
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and the product of the precession and drift wave frequencies 
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is therefore negative in a maximum-J-configuration with a density that increases toward the plasma center, 

' ^da ^ 0- 


( 1 ) 


Typically, trapped-particle instabilities rely on the resonance between these two frequencies, which occurs due to so-called 
“bad” curvature. In quasi-isodynamic configurations, however, we see that trapped particles have average “good” curvature and 
thus exert a stabilizing influence. To demonstrate this mathematically, we proceed from the gyrokinetic system of equations in 
ballooning space. 


^ll^llPa + (W - UJda)ga 
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where 4> is the electrostatic potential, Jq is the zeroth order Bessel function of the first kind, ga = fai + ^f^fao denotes the 
non-adiabatic part of the perturbed distribution function, and the equilibrium distribution function fao is Maxwellian. The ratio 
between the temperature and density gradients is denoted by = d\nTa/d In Ua, and we have written = a;»a[1 + Va— 
3/2)], with = TTiaV^/2Ta. The system of equations is closed by the quasi-neutrality condition. 


= [ gaJod^v. (3) 
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Our argument is based on the energy budget of the instability jstj. We define a quantity 

Pa = eaIm{(w||V||Pa - 0ddaga)4>* Jo} , (4) 

which is the rate of gyrokinetic energy transfer from the electrostatic field to species a. For compactness we have used the 
notation 


{...} 






where A = v]_/v^B and a = '(;||/|'(;|| |. We thus multiply the gyrokinetic equation (|2|i by eaJo(j)*, sum over all species, integrate 
over velocity space and along the entire field line in ballooning space, —oo < I < oo, and take the imaginary part. We note that 
with a complex mode frequency uj = uir + i'y 
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where b = k^Ta/nia^a 
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u < 1. 


Thus we obtain a relation that describes the energy budget of the fluctuations 

o2 
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which is the generalization to an inhomogeneous plasma (in ballooning space) of Eq. (FIO) in @]. The right-hand side represents 
the total energy input from the fluctuations into the various species and must be negative for a growing instability. 
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Now consider a species a with a bounce frequency Wba far above the mode frequency, uj <C uJba, e.g., the electrons in the case 
of ordinary TEMs or both species in the case of the collisionless trapped-particle instability of Rosenbluthj^l and Kadomtsev 
and Pogutse m. We further assume that 0 < ? 7 a < 2/3 (so that a;*a and w// have the same sign for all energies) and that 

1 /2 

uJda has the same sign for all orbits. Thus ordering w ^ u;*o "C k\\ {Ta/ma) we can expand the distribution function, 
ga = 9 a 0 +gal -and obtain 


gaO = 


CaJotj) UJ - 
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and 
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Here we have neglected the passing particles, whose response is a factor a;/A:||t>ra ^ 1 smaller than that of the trapped ones. 
Substituting these results in the expression (|4]i for the energy transfer gives 
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Finally, we consider the limit of this expression when marginal stability is approached, 7 —0+, where we obtain 


Pa = ^ {5(CU - UdapdaicJda - | Jo</>| VaO } • (7) 

J- a 

If and uJda are of opposite signs. Pa > 0 and energy flows from the electric held fluctuations to plasma species a, which 
therefore exerts a stabilizing influence. Consequently, for instabilities with such low frequencies that oj <C oJba for all species, 
we And that > 0, which is in contradiction to Eq. Q at the point of marginal stability, implying the non-existence of a 

marginal stability point and consequently the absence of an instability. The case where the real part of the frequency vanishes, 
ujr — 0 , requires slightly more care, since the resonance then occurs at zero energy and consequently all Pa = 0 , so that an 
instability cannot be ruled out by Eq. 0 . However, from Eq. ®, we obtain for cUr = 0 
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where for small but finite 7 all terms are positive, again in contradiction to Eq. 0; therefore a mode with = 0 at marginal 
stability cannot exist. Hence we conclude that the collisionless trapped-particle mode is absent - i.e., there is no instability with 
frequency far below the ion bounce frequency. This conclusion is an extension of the result by Rosenbluth to an arbitrary number 
of particle species, finite k_\_pa, finite temperature gradients up to pa < 2/3, and finite values of uj/ujda- 

If only the electrons have a bounce frequency that exceeds w but uj ^ ujbi, then we cannot rule out instability by this argument, 
but we can say something about the nature of a possibly occuring mode. We proceed from the gyrokinetic equation 0 and treat 
it as we did to obtain Eq. 0 , only that we do not sum over the species. We then find at the point of marginal stability 


Pa — tU^Im} — ^rQa- 


The distribution function ga app earing in this quadratic form Qa can be obtained from the solution of the gyrokinetic equation 
0 given in Refs. El and ill (correcting for misprints). In the region of velocity space corresponding to trapped particles, 
A > 1 /i3max, where i?niax denotes the maximum field strength along the field line, the solution is 


'^gapl) 


2ea/a0 ^ ^*a 

Ta sin(M(w,Zi,Z2)) 
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where li{X) and hiX) are the bounce points (defined by XB = 1) immediately surrounding I, and where we have written 
lu = max(^, I') and k = min(Z, I') and defined 
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If the growth rate 7 is taken to be positive the solution in the untrapped region is given by 




dt 


dV 


—00 ^ ^ J—00 7II 


4>Jq cos (M(<, I, I')). 


(9) 


The quadratic form Qa can be written as a sum of the contributions from trapped and passing particles of each species 
separately, Qa = Qat + Qap- Substituting the solution (|9]l for passing particles gives 
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where VTa = {“^TalrnaY^^ is the thermal velocity and 
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Therefore, at marginal stability, where cv has an infinitesimal positive imaginary part, we have 
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j=cos,sin 


In ballooning space, there is an infinity of trapping regions along the held line, which are periodic in a tokamak but irregularly 
distributed in a stellarator. When calculating the contribution from the trapped particles to the quadratic form Qa we need to 
sum over all these trapping wells, and then obtain 
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Near marginal stability, an imaginary contribution arises due to the zeros of the sine in the denominator, and when splitting the 
cosines symmetrically we hnd 
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with 


A{x,X,w)=[ dZ(/)Jo cos(M(a;,Zi,0)- 

Jh VI — XB 

We now note that for all species a the forms Qat and Qap have the character of a weighted average over a; of (w — due to 
the positive-dehniteness of the other factors. Thus we can write 

poo 

Pa =uj (a; - a;J^^)PoSa(a:, uj)dx, (14) 

^0 

where PoSa(x, ui) is a positive dehnite function. If we now assume the mode travels in the electron diamagnetic direction, i.e. 
a;w*e > 0, we know from Eq. (|7|i that Pe = 0 due to the lack of resonance. Consequently Eq. (|5]l then implies that Pi = 0 at the 
point of marginal stability. However, from Eq. (fl4l l we obtain Pi > 0 since ww*i < 0, which again implies the non-existence of 
the marginal stability point and the absence of this particular mode. Thus any unstable mode that could arise with uj ^ ujbi must 
propagate in the ion direction at marginal stability and as a consequence from Eq. (|7| draw energy from the ions rather than the 
electrons (Pe > 0 follows from Eq.ljT]), with Eq. (|5]l then implying Pi < 0). There are thus no ordinary TEMs, which tend to 
cause much of the transport observed in tokamaks. We also note that in the usual treatment of the “ubiquitous” mode of Coppi 
and Rewoldt ifoll a;*a • 'ixda < 0 implies stability as well. These conclusions hold as long as 0 < 77 a < 2/3 for all species, and 
collisions may be ignored, but the dissipative TEM could still be unstable. 
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Since this argument is essentially only based on the requirement of quasineutrality and an analysis of the energy budget, it 
is independent of all geometric details of the magnetic field except the condition that the bounce-averaged curvature should be 
favorable, dJ/dip < 0, for all orbits. This requirement can also be satished in other omnigenous conhgurations lll4i Il5|] . In 
a tokamak, for example, it is achieved if the pressure gradient is steep enough to cause drift reversal of all trapped particles 
Olhll . but in practice such a steep pressure gradient necessitates taking account of electromagnetic effects. However, if MHD 
ballooning modes are stabilized by negative magnetic shear (according to the tokamak definition), it is conceivable that the 
stabilization of trapped-electron modes could help explain the transport reduction observed in internal transport barriers. 

It is, of course, difficult to achieve exact quasi-isodynamicity, but one expects that the drive for trapped-particle modes should 
become weak if most orbits satisfy UJ^,e ■ oJde < 0. One would expect that even approximately quasi-isodynamic stellarators 
should have relatively small trapped-particle instability growth rates, particularly if central fueling is accessible through pellet 
injection so that a stabilizing density gradient can be achieved. Finally, it should be mentioned that inverting the density gradient 
in a tokamak, so as to reverse the sign of a;*e, has long been known to make the collisionless trapped-electron mode less unstable 
I ITlllSl], because there are then fewer electrons with • uJ^ > 0. However, this stabilization is incomplete since in a typical 
tokamak there are always electrons with both signs of 

In summary, whereas in tokamaks most of the transport in the core tends to be driven by ITG and ETG modes, and by modes 
driven unstable by trapped electrons, the latter are stable in quasi-isodynamic stellarators in the electrostatic and collisionless 
limit, if 0 < 77 q < 2/3, and so are also all such instabilities with frequencies below the ion bounce frequency. 
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